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Permanent focal brain damage can have critical effects on the function of nearby as 
well as remote brain regions. However, the effects of transient disturbances on global 
brain function are largely unknown. Our goal was to develop an experimental in vivo 
model to map the impact of transient functional brain impairment on large-scale neural 
networks in the absence of structural damage. We describe a new rat model of 
transient functional hemispheric disruption using unilateral focal anesthesia by intracarotid 
pentobarbital injection. The brain's functional status was assessed with resting-state fMRI 
(rs-fMRI) and electroencephalography (EEG). We performed network analysis to identify 
and quantify highly connected network hubs, i.e., "rich-club organization," in pre- and 
postbarbital functional networks. Perfusion MRI data demonstrated that the catheterized 
carotid artery predominantly supplied the ipsilateral hemisphere, allowing for selective 
hemispheric brain silencing. The prebarbital baseline network displayed strong functional 
connectivity (FC) within and between hemispheres. Following pentobarbital injection, the 
disrupted hemisphere revealed increased intrahemispheric FC with concomitant decrease 
of interhemispheric connectivity. The bilateral functional network was characterized 
by a strong positive rich-club effect, which was not affected by ipsilateral disruption. 
Nevertheless, the rich-club value was significantly decreased in the ipsilateral hemisphere 
and to a lesser extent contralateral^. Loss of interhemispheric EEG synchronization 
supported the rs-fMRI findings. Our data support the concept that densely connected 
rich-club regions play a central role in global brain communication, and show that network 
hub configurations can be significantly affected by focal temporary functional hemispheric 
disruption without structural neuronal damage. Further studies with this rat model will 
provide essential additional insights into network reorganization patterns in response to 
transient functional brain disruption. 
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INTRODUCTION 

Focal brain damage, such as after ischemic stroke or in epilepsy, 
not only affects local neuronal networks, but can also have signifi- 
cant effects on the function of remote brain regions. For instance, 
constrained behavioral output (Carter et al., 2010) and cogni- 
tive comorbidities (Tatemichi et al., 1994; Motamedi and Meador, 
2003) after stroke and epilepsy have partly been attributed to 
functional disturbances of distant brain regions that were not 
primary affected. In addition, spontaneous partial recovery of 
sensorimotor and behavior output in the weeks and months fol- 
lowing stroke may be explained by functional plasticity of remote 
and unaffected brain regions (Langhorne et al, 2009; Stinear, 
2010). As brain functions are believed to emerge from interacting 
regions within large-scale neural networks, the brain's intrinsic 
capacity to compensate for focal injury is linked to reorganiza- 
tion of surviving functional networks at a global level (Grefkes 
and Fink, 2011). 

A sophisticated set of measures, using the concept of graph 
theory, has been used to describe alterations in functional 



network organization following focal brain damage (Rubinov and 
Sporns, 2010), with promising opportunities for applications in 
clinical settings (Fornito et al, 2013). Obviously, changes in func- 
tional networks are strongly dependent on the structural status of 
neuronal connections that link brain regions. Recent studies have 
reported a high similarity between functional network connectiv- 
ity and underlying structural neuronal connections (Achard et al., 
2008; Honey et al., 2009). However, it is largely unknown how 
functional networks reorganize in response to transient focal dis- 
turbances without structural damage. Characterizing the impact 
of transient localized functional brain disruptions on large-scale 
neural networks can provide important insights into a network's 
resilience, which might lead to improved understanding of under- 
lying causes of comorbidities in diseases with localized brain 
damage, such as focal epilepsy and ischemic stroke. 

Global neural networks are known to have a hierarchical struc- 
ture with the existence of a number of well-connected and highly 
central hub regions (Hagmann et al, 2008; Bullmore and Sporns, 
2012) that form a so-called "rich-club," which reflects a network's 
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capacity to segregate and integrate multisensory information in 
a united manner (Van Den Heuvel and Sporns, 2011; Harriger 
et al, 2012). These well-connected regions play a key role in global 
signal processing between different parts of the network, and 
functional disturbance of hub regions can have critical impact on 
overall network functioning. In this study we aimed to map the 
impact of transient functional brain disruption on global neural 
networks in the absence of structural damage in a newly devel- 
oped reproducible rat model. To that aim, we characterized the 
brain's functional network organization with resting-state fMRI 
(rs-fMRI) after unilateral pentobarbital-induced brain "silenc- 
ing." We specifically focused on effects on well-connected hub 
regions using a relatively new framework that provides quantita- 
tive measures of the organization of highly connected rich-clubs 
(Van Den Heuvel and Sporns, 2011; Harriger et al, 2012; Collin 
et al., 2013). So far it has remained unknown whether rodent 
brain networks contain rich-club connections and if so, how 
these connections should be interpreted. Our hypothesis was 
that (i) rat brains have functional rich-club connections, and (ii) 
that transient hemispheric functional disruption affects the rich- 
club organization in both the disrupted and unaffected network 
regions. 

METHODS 
RAT MODEL 

We developed a rat model of transient functional hemispheric 
disruption using unilateral focal anesthesia, similar to the intra- 
carotid sodium amobarbital procedure for the "Wada test" in 
humans (Wada and Rasmussen, 2007). This procedure was cho- 
sen to temporally disrupt intra- and interhemispheric functional 
connectivity (FC) of the ipsilateral hemisphere. The protocol was 
approved by the Ethical Committee on Animal Experiments of 
the University Medical Center Utrecht and Utrecht University. 
Experiments were carried out in accordance with the guidelines 
of the European Communities Council Directive. 

Three healthy adult Sprague-Dawley rats (Charles River 
Laboratories International Inc., MA, USA) were mechanically 
ventilated with 2% isoflurane in air: 02(7:1; rate 52-59min~ 1 ). 
To catheterize the right carotid artery, the right common, and 
external carotid arteries were exposed by a ventral incision in 
the neck. A small catheter was advanced into the external carotid 
artery in the direction of the common carotid artery until the 
blunt catheter tip reached the bifurcation of the internal and 
external carotid arteries. Blood flow from the common into the 
internal carotid artery was ensured to be maintained (that is, the 
catheter tip remained in the ligated external common artery). The 
catheter was kept in place with a 5-0 suture and instant glue. 
The incision was closed and animals received 5 mL glucose-saline 
solution to compensate for loss of water and minerals. 

Directly after surgery, rats were placed in an MR-compatible 
stereotactic holder and immobilized with earplugs and a tooth- 
holder. Blood oxygen saturation and heart rate were continuously 
monitored with a pulse-oximeter with the probe positioned on a 
hindpaw. In addition, expired CO2 was continuously monitored 
with a capnograph, and body temperature was maintained at 37.0 
± 0.5° C using a feedback-controlled heating pad. End-tidal CO2 
levels were kept within the normal range, equivalent to arterial 



pCC>2 levels between 35 and 45 mmHg (calculated from previous 
calibration measurements), by adjusting ventilation volume and 
rate when necessary. 

MRI 

MRI experiments were done at 4.7 T with a 125 mm internal 
diameter gradient-coil insert. A Helmholtz volume coil and an 
inductively coupled surface coil were used for signal excitation 
and detection, respectively. 

Repetitive blood oxygenation level-dependent rs-fMRI scans 
were acquired using a gradient echo single shot echo-planar 
imaging (EPI) sequence [repetition time (TR)/echo time (TE) = 
500/19 ms; 7 1.5-mm coronal slices; field-of-view = 32 x 32 mm; 
in-plane voxel resolution, 0.5 x 0.5 mm; pulse angle = 35°; tem- 
poral scan resolution = 500 ms; total scan time = 15 min] . Exactly 
10 min prior to rs-fMRI acquisition, end-tidal isoflurane con- 
centration was reduced to 1.0%, which was maintained during 
rs-fMRI. 

After 7.5 min rs-fMRI acquisition, 175 (xL of 2mg/kg pen- 
tobarbital was slowly injected intracarotidly. This concentration 
was based on an electroencephalography (EEC) titration/dose- 
response experiment, which showed that suppression of EEC 
activity occurred at this threshold level at 1% isoflurane anes- 
thesia (data not shown). EEC activity returned to baseline within 
15 min. 

In one animal, perfusion MRI was conducted (TR/TE = 
100/25 ms) following rs-fMRI, using gradient echo EPI (TR/TE = 
100/25 ms; 11.5-mm coronal slice; field-of-view = 32 x 32 mm; 
in-plane voxel resolution, 0.5 x 0.5mm; pulse angle = 90°) in 
combination with an bolus injection of 50 (xL of l.Ommol/mL 
gadobutrol through the intracarotid catheter. Serial images were 
inspected to verify unilateral bolus arrival. 

EEG 

To validate functional MRI measurements, EEG was recorded 
immediately following rs-fMRI acquisitions, at 1.0% isoflurane 
anesthesia outside the MR scanner. To that aim a small medial 
incision was made in the skin covering the skull and the peri- 
cranium. A hole in the skull above the left and right primary 
somatosensory cortices was made with a 300 \im micro drill. 
Next, epidural EEG was recorded from the bilateral somatosen- 
sory cortices with a reference electrode above the cerebellum. 
EEG activity was recorded for 15 min in two animals. After 
7.5 min of EEG measurements, 175 /xL of 2mg/kg pentobar- 
bital was slowly injected intracarotidly. The time between this 
injection and the first injection in the MR scanner was > 1 h. 
EEG measurements were conducted using a homebuilt multi- 
channel amplifier, and band-pass filtered between 0.1 Hz and 
250 Hz, using a National Instruments™ NI USB-6211 DAQ with 
a sampling rate of 1000 Hz per channel. 

IMAGE PROCESSING AND ANALYSIS 

Rs-fMRI images were motion-corrected, spatially smoothed, 
band-pass filtered between 0.01 and 0.1 Hz, and nonlinearly 
matched with the Paxinos and Watson rat brain atlas (Paxinos 
and Watson, 2006) for FC analysis. The first (prebarbital) and last 
(postbarbital) 500 images were used to extract time series for 32 
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regions-of-interest (ROIs) in each hemisphere, and to calculate 
pair-wise FCs between all ROIs from the correlation coefficient r. 
ROIs were carefully selected to basically include all bilateral neo- 
cortical and subcortical gray matter atlas regions (see Figure 2). 
The cerebellum was excluded from analysis. All ROIs were large 
enough to contain multiple voxels at the applied fMRI resolution, 
and were functionally connected to at least one other region. 

ESTIMATION OF WEIGHTED RICH-CLUB EFFECT 

Pre- and postbarbital FCs were averaged across animals to con- 
struct two 64 x 64 node networks (each node is a cortical or 
subcortical ROI in one of the hemispheres). From these net- 
works ipsilateral and contralateral 32 x 32 node sub-networks 
were extracted. Negative FCs were excluded to allow calculations 
according to the weighted rich-club methodology (Opsahl et al., 
2008). 

The weighted rich-club effect characterizes the tendency of 
network hubs to engage in stronger or weaker FCs among them 
than expected by chance. The calculation involves multiple steps. 
First, all network nodes are ranked in terms of prominence. We 
used node strength s as a measure of prominence. Strength is 
defined as the sum of FCs attached to the connections originating 
from a network node (Barrat et al., 2004). For increasing values 
of s, we selected the group (the club) of all nodes whose strengths 
were larger than s. In this step we obtain a series of increasingly 
selective clubs. For each of these clubs we counted the number 
E >s of connections connecting the club members, and measured 
the sum W >5 of the FCs attached to these connections. 

In the next step we determined the ratio 4> w (s) between W >s 
and the sum of the FCs attached to the E >s strongest connections 
within the network. Mathematically: 



where w\ ank > vf[™* with / = 1, 2, . . . , E the ranked FCs on the 
connections of the network. E is the total number of connections. 
With this formulation we determined the fractions of FCs shared 
by rich nodes compared with the total amount they could share 
if they were connected through the strongest connections of the 
network. 

A high value of (p w (s) is not sufficient to account for the 
weighted rich-club effect as higher- strength nodes in a network 
have a higher probability of sharing connections with each other 
just by chance. In other words, networks where connections are 
randomly established could have a non-zero value of (p w {s). To 
control for this effect we normalized 4> w (s) relative to the aver- 
age weighted rich-club effect, 0 r " n d om ( s )> found in 1000 random 
networks generated using a reshuffle procedure. This procedure 
reshuffles the FCs globally in the networks, while keeping the 
topology intact. 

The weighted rich-club effect p w (s) at strength s is the ratio: 
p w {s) J W{s) , When p w {s) is > 1.0 the functional network 

^random ' S ' 

has a positive weighted rich-club effect, with nodes in the rich- 
club concentrating most of their connections toward other club 
members relative to random networks. 



EEG ANALYSIS 

EEC signals were analyzed in consecutive epochs of 30 s. For each 
epoch the correlation coefficients between the temporal fluctua- 
tions of the EEG signals from left and right somatosensory cor- 
tices were calculated for the delta (1-5 Hz), theta (5-8 Hz), alpha 
(8-15 Hz), beta (15-30 Hz), gamma (30-250 Hz), and broad 
(1-255 Hz) frequency bands. 

RESULTS 

Perfusion MRI data demonstrated that the catheterized carotid 
artery predominantly supplied the ipsilateral hemisphere allow- 
ing for selective hemispheric brain silencing with pentobarbital 
(Figure 1). Perfusion was strongest in the cortical areas, whereas 
levels in subcortical regions were lower. We also detected some 
contrast agent arrival in the anterior contralateral hemisphere, 
which can be explained by delivery via the circle of Willis. 

NETWORK PROMINENCE 

The prominence of all network nodes, measured as the strength, 
is plotted in Figure 2 for the prebarbital (baseline) and post- 
barbital whole-brain functional networks, as calculated from 
the rs-fMRI experiments. The ten most prominent prebarbital 
ROIs included mostly somatosensory regions, namely, the ipsi- 
and contralateral dysgranular zone of the primary somatosen- 
sory cortex; the contralateral forelimb region, upper lip region 
and barrel field of the primary somatosensory cortex; and 
the ipsilateral jaw region of the primary somatosensory cor- 
tex. The agranular and granular insular cortex and the facial 
nuclei were also prominent network nodes before pentobarbital 
injection. 

Following unilateral pentobarbital injection, a subgroup of 
these ROIs lost prominence (Figure 2). This was most pro- 
nounced for the contralateral barrel field of the primary 
somatosensory cortex and the agranular insular cortex. In con- 
trast, the prominence of the contralateral primary motor cortex 
increased toward the highest value in the postbarbital network. 
The ROIs presenting with the largest changes in prominence were 
not confined to a single anatomical region. 

FUNCTIONAL CONNECTIVITY 

The average baseline functional network, i.e., before pento- 
barbital injection, displayed strong connectivity within and 
between hemispheres. After pentobarbital injection the ipsilat- 
eral hemisphere revealed increased intrahemispheric connections 
with concomitant decrease of interhemispheric FC (Figure 3 A). 
Increased intrahemispheric FC was found in both hemispheres. 
Several ipsilateral regions retained strong network strengths with 
other regions, including the facial nucleus, prelimic cortex, cin- 
gulate cortex area 1, and the granular insular cortex. 

WEIGHTED RICH-CLUB EFFECT 

The bilateral functional network was characterized by a strong 
weighted rich-club effect, with p w (s) > 1.0 over a range of 
prominence thresholds (Figure 3B). This rich-club effect was not 
affected by ipsilateral FC disruption. In contrast to the whole 
brain network, the p w (s) values were significantly decreased 
in the ipsilateral sub-network following unilateral pentobarbital 
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FIGURE 1 | Perfusion MRI frames, demonstrating clear unilateral delivery top left to bottom right). The frames are in coronal plane at 1 mm from bregma, 

of the contrast agent in the hemisphere ipsilateral to the side of Presence of contrast agent is displayed using a color scale: dark for low 

cannulation. The frames were acquired at 100 ms intervals from the start of contrast concentrations and bright for high contrast concentration (in arbitrary 

bolus administration until 1.5 s later (the temporal frame sequence runs from units). The green rectangles mark the ipsi- and contralateral hemispheres. 



injection (Figure 3C). This pattern of decreased p w (s) values was 
to a lesser extent also seen in the contralateral sub-network. 

DATA VARIATION 

The variation between individual functional networks is shown 
in Figure 4. The standard deviation of the FC values slightly 
increased post pentobarbital injection as shown from the density 
plots. The standard deviation range was 0-0.3 prior to injection 
and 0-0.4 post injection. 

ELECTROENCEPHALOGRAPHY 

The EEG measurements, acquired in the same animals directly 
after imaging, supported our MRI -based observation of decreased 
interhemispheric resting-state FC. Interhemispheric EEG syn- 
chronization decreased after unilateral pentobarbital injection in 
the alpha, theta, and broadband frequencies (Figure 5). The syn- 
chronization in the other frequency bands increased temporarily 
(gamma) or remained unchanged (beta). 

DISCUSSION 

Our pilot study shows that transient functional global network 
disruptions in response to unilateral and reversible functional 
brain disturbance can be effectively studied using a straightfor- 
ward rodent model in combination with rs-fMRI. Our rat data 
support the concept that densely connected rich-club regions 
play a central role in global brain communication, as has been 
previously shown in humans. In addition, we found that global 
network hub configurations can be significantly affected by tem- 
porary functional disruption of a hemisphere. 



A rodent model that allows selective and reversible functional 
hemispheric silencing has not been described before. Our results 
show that this model can be useful to assess effects of (tran- 
sient) brain impairment on functional networks. The concept 
of this model was based on recent FC and network data pre- 
sented by Douw et al. (2009, 2010). They showed that FC within 
and between hemispheres changed in persons with epilepsy dur- 
ing an intra-arterial amobarbital injection — the WADA test — 
as measured with EEG (Douw et al, 2009). They concluded 
that the effects of unilateral anesthesia on FC in the contralat- 
eral hemisphere should not be ignored. Direct comparison of 
low-frequency rs-fMRI connectivities with the different EEG fre- 
quency bands is difficult, and more data are required to elucidate 
differences and similarities between rs-fMRI and EEG patterns. 
Nevertheless, the observed decrease in interhemispheric EEG FC 
in the delta and theta bands by Douw and colleagues corre- 
sponds with our interhemispheric EEG results (Figure 5). Over 
the ipsilateral hemisphere the delta and theta band synchro- 
nization increased after amobarbital injection in humans, which 
also agrees with our findings of increased ipsilateral intrahemi- 
spheric FC in rat brain after pentobarbital injection. Apparently, 
unilateral functional brain disruption by amo- or pentobarbi- 
tal injections is confined to subverted functional outcome (that 
is, speech and language comprehension in humans) rather than 
absolute inhibition of low-frequency FC. This idea is further 
supported by recent FC measurements at an anesthetic depth 
characterized by unresponsiveness, where reductions in con- 
nectivity were shown to be partial and not complete (Hudetz, 
2012). 
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FIGURE 2 | Pre- and postbarbital individual network node strengths 
plotted as bar graphs. The 64 bilateral nodes are sorted from high to 
low strength. The 10 network nodes with the highest strength prior to 
pentobarbital injection are given in blue. The node with the highest 
strength in the postbarbital functional network, the contralateral primary 
motor cortex, is given in orange. The following abbreviations are used, 
contra, left hemisphere, contralateral to cannulation, ipsi, right 
hemisphere, ipsilateral to cannulation; AID, agranular insular cortex dorsal 
part; AIP, agranular insular cortex posterior part; Aul, primary auditory 
cortex; AuD, secondary auditory cortex dorsal area; AuV, secondary 
auditory cortex ventral area; AIV, agranular insular cortex ventral part; Cgl, 



cingulate cortex area 1; Cg2, cingulate cortex area 2; CPu, 
caudate-putamen complex; Dl, facial nucleus; Fr3, frontal cortex area 3; 
Gl, granular insular cortex; M1, primary motor cortex; M2, secondary 
motor cortex; PrL, prelimbic cortex; PtPD, parietal cortex posterior area 
dorsal part; PtPR, parietal cortex posterior area rostral part; Rt, 
intermediate reticular nucleus; S1BF, primary somatosensory cortex barrel 
field; S1DZ, primary somatosensory cortex dysgranular zone; S1FL, 
primary somatosensory cortex forelimb region; S1HL, primary 
somatosensory cortex hindlimb region; S1J, primary somatosensory 
cortex jaw region; S1UI_p, primary somatosensory cortex upper lip region; 
S2, secondary somatosensory cortex; TeA, temporal association cortex. 



Douw et al. also assessed EEG data acquired during WADA- 
test using whole-brain functional network analysis (Douw et al., 
2010). Their results revealed an increased random network orga- 
nization after unilateral functional disruption with amobarbital 
in people with epilepsy (Douw et al., 2010). However, the network 
topologies within the injected and non-injected hemispheres were 
not studied separately in this group of patients. In our study 
we found no difference in rich-club characteristics between pre- 
and postbarbital whole-brain networks. Differences in weighted 



rich-club effects were only found in sub-network analyses of indi- 
vidual hemispheres (Figure 3). The asymmetric presentation of 
prebarbital hub regions (that is, mostly contralateral sensorimo- 
tor areas) is in line with studies in healthy rats that reported 
right-sided population bias and functional lateralization (Glick 
and Ross, 1981; Alonso et al, 1991). 

A positive weighted rich-club effect characterizes the tendency 
of network hubs to engage in stronger functional connectivities 
among them than expected by chance. We found significantly 
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FIGURE 3 | Pre- and postbarbital (A) whole-brain functional connectivity 
(FC) networks, displayed as matrices. Each row and column stands for a 
region-of-interest in the hemisphere ipsi- or contralateral to the side of 
cannulation. Matrix elements are colored according to FCs values, measured 
as the correlation coefficient r between regions-of-interest. Red indicates no 
FC; yellow indicates high FC. Interhemispheric FCs are delineated with blue 
squares; intrahemispheric FCs with green squares. Weighted rich-club 
effects (y-axis) are plotted against increasing levels of strength (x-axis) for the 



entire brain network (B), and the contra- and ipsilateral (C) hemispheric 
networks. The maximum strength threshold depends on the size of the 
network. The strength thresholds were increased until networks started to 
disintegrate. This explains the different scaling of the x-axis in the 
hemispheric networks as compared to the entire brain network. Significant 
differences in weighted rich-club effects between prebarbital (blue line) and 
postbarbital (red line) networks are indicated as * [permutation analysis; 
p < 0.01 considered significant]. 



decreased p w (s) values in the ipsi- and contralateral sub-networks 
after unilateral pentobarbital injection. Similar positive rich-club 
organizations have been described at a neuronal level in the 
nematode worm Caenorhabditis elegans (Towlson et al., 2013) and 
at a larger scale in human neocortical networks (Van Den Heuvel 
and Sporns, 2011). Our findings in rodent brain networks, con- 
solidate the hypothesis that connections according a rich-club 
layout may be a general and scale-invariant principle of brain 
network organization (Towlson et al, 2013). Nonetheless, further 
validation of this hypothesis is essential, given the small sample 
size used in this study. 

The rich-club connections in the structural human brain 
network are considered to be the backbone for global brain 
communication (Van Den Heuvel et al., 2012). The majority 
of these connections are long-distance neural pathways, such 
that the rich-club phenomenon is costly in terms of consumed 
energy and space. However, the trade-off for higher cost is higher 



performance in terms of network integration; up to 70% of all 
network's shortest paths travel through rich-club connections. 

Whether structural network principles also hold for functional 
rich-club networks remains to be determined, but we expect 
the functional rich-club organization to follow the underlying 
structural organization to a large extent. On the other hand, 
functional rich-club connections may be much more flexible in 
acutely responding to environmental changes, such as induced 
by anesthesia. This is, for instance, underscored by the differ- 
ences in our whole-brain and sub-network data (Figures 3B,C). 
Unilateral brain silencing did not affect the rich-club effect at a 
whole-brain level; only the club members changed dramatically 
(Figure 2). On the other hand, the hemispheric sub-networks 
did reveal change in rich-club connectivity after hemispheric 
silencing. Apparently the level of network definition is criti- 
cal for detecting changes in functional rich-club organization 
after functional brain lesioning. A better understanding of this 
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FIGURE 4 | Pre- and postbarbital (A) variation within the whole-brain 
functional connectivity matrices as shown in Figure 3. Each row and 
column represent a region-of-interest in the hemisphere ipsi- or 
contralateral to the side of cannulation. Matrix elements are colored 
according to the standard deviation (SD) of the individual functional 
connectivity values, measured as the correlation coefficient r between 



regions-of-interest. Red indicates SD values close to 0; yellow indicates 
SD values close to 0.4. As shown in the density distribution of all SD 
values (B), the variation between individual networks is slightly smaller 
for the prebarbital data (most SD's between 0 and 0.3) as compared to 
the variation in the postbarbital data (with most SD's between 0 
and 0.4). 



phenomenon may require the acquisition of structural diffu- 
sion tensor imaging data alongside rs-fMRI in future stud- 
ies. This will allow comparison of functional networks and 
their rich-club organization with data from the intact structural 
networks. 

Limitations in this proof-of-principle methods study should 
be acknowledged. First, experiments were conducted at 1.0% 
isoflurane anesthesia. We do not know how this baseline 
anesthesia may have interacted with the pentobarbital-induced 
brain silencing. We expect this effect to be relatively small 
as cerebral glucose utilization, cerebral blood flow coupling 
(Maekawa et al., 1986; Lenz et al., 1998) and FC, includ- 
ing the long distance connections (Liang et al, 2012), are 
largely preserved at 1.0% isoflurane anesthesia (Wang et al., 
2011). We also know that the neurovascular coupling in rat 
brain is similar across distinct types of anesthetics including 
pentobarbital and isoflurane (Franceschini et al., 2010), how- 
ever further validation studies, including robust dose-response 
experiments, are needed. Rs-fMRI acquisition may be ide- 
ally executed in awake rats, however, this is complicated due 
to motion artifacts related to post-surgical pain, stress and 
discomfort. 

Second, anesthesia-based functional brain disruption by uni- 
lateral injection through the carotid artery may not be exclusively 
unilateral due to presence of the circle of Willis and recirculation 



effects. Nonetheless, the small injection volume used in our study 
limited this effect, which was confirmed by the strongly lateralized 
results for hemispheric FCs. 

Third, the likelihood that a statistically significant result 
reflects a true effect is reduced in studies with a small sam- 
ple size (Button et al., 2013). The large changes in functional 
network connectivity and organization that we observed in 
this rat model alleviate this problem. Nonetheless, firm con- 
clusions on directionality of network changes and size of effect 
require additional studies with larger groups of animals. This 
is emphasized by the variation that we found between FC 
correlation coefficients obtained from the individual networks 
(Figure 4). This variation is however comparable to the vari- 
ation found between individual interhemispheric FC in rats 
recovering from a unilateral stroke, as studied with rs-fMRI 
with a similar experimental setup and design (Van Meer et al., 
2012). 

In conclusion, our pilot study suggests that changes in rich- 
club organization following functional network disruption might 
affect the capacity to efficiently segregate and integrate multisen- 
sory information. These changes are unlikely to depend on the 
presence of structural damage. We speculate that functional rich- 
club alterations may reflect an underlying mechanism of neu- 
rological impairments in patients with no or limited structural 
brain damage, such as after transient ischemic attack. 
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FIGURE 5 | The interhemispheric somatosensory functional EEG 
connectivities (measured as correlation coefficients (r)) for consecutive 
epochs, each 30 s, plotted over time (1 5 min) for the distinct frequency 
bands (A). The corresponding moment of pentobarbital injection and 
periods of resting-state fMRI used to construct functional networks are 
shown below (B). Note that functional imaging was acquired prior to EEG, 
but was identical in approach and timing of pentobarbital injection. 
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